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Abstract 

Background: EPR-based distance measurements between spin labels in proteins have become a valuable tool in 
structural biology. The direct translation of the experimental distances into structural information is however often 
impaired by the intrinsic flexibility of the spin labelled side chains. Different algorithms exist that predict the 
approximate conformation of the spin label either by using pre-computed rotamer libraries of the labelled side 
chain (rotamer approach) or by simply determining its accessible volume (accessible volume approach). Surprisingly, 
comparisons with many experimental distances have shown that both approaches deliver the same distance 
prediction accuracy of about 3 A. 

Results: Here, instead of comparing predicted and experimental distances, we test the ability of both approaches 
to predict the actual conformations of spin labels found in a new high-resolution crystal structure of spin labelled 
azurin (T21R1). Inside the crystal, the label is found in two very different environments which serve as a challenging 
test for the in silico approaches. 

Conclusions: Our results illustrate why simple and more sophisticated programs lead to the same prediciton error. 
Thus, a more precise treatment of the complete environment of the label and also its interactions with the 
environment will be needed to increase the accuracy of in silico spin labelling algorithms. 



Background 

The structural characterization of proteins by pulsed 
EPR methods such as pulsed electron electron double 
resonance (PELDOR, also known as DEER) has become 
increasingly popular in recent years [1-3]. A prerequisite 
for the PELDOR experiment is the presence of at least 
two paramagnetic centres in the protein. Since most 
proteins are diamagnetic, spin labels such as MTSSL [4] 
are routinely attached to the surface of proteins via site- 
directed spin labelling techniques [5-7]. The PELDOR 
experiment measures distances between such spin labels. 
The obtained distance information can then be used to ana- 
lyse conformational changes e.g. of membrane proteins 
[8-10] or to reconstruct macromolecular complexes [11,12], 
Although spin labels, such as MTSSL, are small compared 
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to FRET labels (MTSSL has roughly the size of an arginine 
residue), they still act as a flexible link between the protein 
and the spin-centre itself (usually the nitroxide group of 
MTSSL). Thus, the problem arises, that on the one hand, 
the PELDOR experiment delivers an accurate distance 
between the two spin centres, but on the other hand, the 
exact position of the distance vector with respect to the 
protein is unknown. Many studies have investigated this 
problem. For example, the crystallographic analysis of 
multiple spin labelled T4-lysozyme mutants provided in- 
sights into the interactions of spin labels with proteins 
and preferred rotameric states of the labels [13-16]. The 
problem also led to the development of in silico spin label- 
ling programs, which aim to predict the conformation of 
the spin label in the local environment of the attachment 
site [17-21]. These programs attach a model of the spin 
label to a macromolecular structure and generate an en- 
semble of possible rotamers of the label. The programs 
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either allow all possible rotamers (e.g. for MTSSL: the five 
dihedral angles (xl- X^) are randomly set, [20]) or the 
angles are set based on calculated rotamer libraries and/or 
crystal structures [18,19,22]. With each method, distances 
between the generated ensembles can be measured and 
compared to experimental distances. Surprisingly, exten- 
sive benchmarks have revealed that independent of the 
approach, the average error of the predicted distances is 
around 3 A, [20,22,23]. This indicates that current in silico 
spin labelling programs do not accurately enough model 
the complex interactions between a spin-label and its 
molecular environment. It is therefore important to gain 
further experimental insight into label-environment inter- 
actions. This will help to further extend available rotamer 
libraries with experimental rotamer data and provide clues 
concerning possible improvements of in silico spin label- 
ling algorithms. 

We present here the high-resolution (1.5 A) structure 
of azurin from Pseudomonas aeruginosa PAOl spin 
labelled with MTSSL at position T21. Due to the crystal- 
packing environment in the triclinic crystals, the label is 
observed in two different but fully ordered states. The 
excellent quality of the electron density allows precise 
measurements of all five dihedral angles in both con- 
formations. Whereas one of the conformations fits to 
calculated rotamer libraries [18] the second conformation 
does not, illustrating how the environment of the label 
can force it to adopt a conformation which is less 
favourable for the free label. The implications for in silico 
spin labelling are discussed. Further, we present a newly 
developed, affinity chromatography based purification and 
spin labelling protocol for azurin. 

Results 

Protein production, spin labelling, crystallisation and 
structure solution 

The classic purification protocol for azurin describes the 
isolation of recombinant azurin from the periplasmic 
fraction of E. coli [24]. Although this reportedly works 
well, we aimed to use affinity chromatography methods 
to streamline the production of azurin and to be able to 
use "on-column" spin labelling. To achieve this, the cod- 
ing sequence for the mature azurin protein (excluding 
the coding sequence for the signal sequence) was cloned 
into the pHisGSTTEV vector. The T21C mutant was 
introduced by PCR techniques [25] . The protein was then 
overexpressed in E. coli, bound to the Ni 2+ affinity resin, 
washed and spin labelled during the elution step. The 
GST tag was cleaved off with TEV protease and removed 
by gel filtration. Using this method, ~1 mg of pure, spin 
labelled azurin was obtained from a 1 1 bacterial culture. 
Note that this yield is lower than the yield reported for the 
conventional protocol [24]. However, in our hands both 
methods yielded very similar amounts of pure azurin. The 



spin labelled protein (10 mg/ml) was crystallized at room 
temperature using a commercial crystallisation screen 
(JCSG+) and several initial hit conditions were obtained. 
The components of these conditions were combined in a 
stochastic optimisation screen. After several days, plate 
shaped, blue azurin crystals (-100 um) appeared. A 200° 
diffraction dataset was collected at BESSYII (BL14.1). The 
data could only be processed and scaled in space group 
PI. The structure was then solved by molecular replace- 
ment using the structure of zinc bound azurin as search 
model (PDB-ID: 1E67 [26]) and refined until R/R free con- 
verged at 18.8/21.4. Data collection and refinement statis- 
tics are summarized in Table 1. The refined structure was 
analysed with LABELIT to confirm that no higher crystal- 
lography symmetry than PI is present [27]. 

Overall structure of azurin T21R1 

The asymmetric unit of spin labelled azurin T21R1 
contains four independent monomers of azurin (I-IV) that 
are related by a two-fold NCS axis (Figure 1). The electron 
density for all four monomers is of excellent quality and a 
complete model with very good stereochemistry could be 

Table 1 Data collection and refinement statistics 



Data collection 

PDB-ID 
Space group 
Unit cell (A, °) 

Matthews coefficient (A 3 Da" 1 ) 
Solvent content (%) 
Molecules per ASU 
Wavelength (nm) 
Beamline 

Resolution range (A) 
Total observations 
Unique reflections 
Completeness 
Multiplicity 

Emerge 
Rmeas 

l/o(l) 

Refinement 
Resolution range (A) 

R/Rfree (%) 

RMSD bonds (A)/angles (°) 

Ramachandran plot, 

(% fa vo u red/a 1 1 o wed/d i sa 1 1 o wed) 

MOLPROBIPr' score/clash score 



Azurin T21R! 

4BWW 
P1 

37.0, 53.7, 73.2, 74.4, 89.3, 83.4 
2.41 
49.0 
4 

0.918409 
BL14.1, BESSYII 
51.3-1.48 
155295 
82307 
91.8 (90.9) 

1.9 (1.9) 
0.07 (0.47) 
0.1 1 (0.88) 
0.07 (0.60) 
5.5 (1.1) 

48.1-1.48 
18.8/21.4 
0.007/1.143 
98.4/1 .6/0 

1 .06/2.77 



Values in parenthesis refer to the shell of highest resolution. 
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Figure 1 Overall structure of triclinic azurin T21R1 crystals. A) The asymmetric unit. The four azurin chains (l-IV) are depicted as cartoon 
models; the R1 side chains are shown as sticks. Brown spheres mark the positions of the copper ions. The two-fold NCS symmetry axis is indicated 
and the monomers are coloured to emphasize the NCS relationship. B) Superposition of all four monomers (colouring is identical to panel A). 
C) Superposition of monomer I with a previously published structure of azurin (magenta, PDB-ID: 1E67, [26]). 



build (Molprobity score/clash score 1.06/2.77, Table 1). All 
four molecules have nearly identical structures amongst 
each other (RMSD = 0.3 A for C a atoms, Figure IB) and to 
other azurin structures (e.g. RMSD = 0.31 A for C a atoms 
to 1E67 structure, Figure 1C). 



glycerol or ethylene glycole (up to 50%) are also used to 
cryo-protect PELDOR samples. Thus, it might well be, 
that in other cases, these molecules influence the dynam- 
ics of the spin label, either by direct interactions or indir- 
ectly by restricting its accessible volume. 



Structure of the R1 side chain 

After the molecular replacement step, positive difference 
electron density for the label was observed at the ex- 
pected positions (T21). Following initial refinement with 
omitted Rl side chain, models of the label could be un- 
ambiguously placed into the electron density in all four 
monomers, resulting in two distinct conformations of 
the label (Figures IB, 2A, B). In monomers I and III, the 
label is involved in vdW interactions with residues T126 
and K128 of the same monomer and vdW interactions 
with residues A119 and L120 of a neighboring azurin 
monomer (Figure 2C). An analysis with the PISA server 
[28] revealed an accessible surface area of 228 A 2 for the 
Rl sidechain at this position, 49 A 2 of which are buried. 
Henceforth, the Rl conformation found in monomers I 
and III will be referred to as Rl-I/III. In monomers II 
and IV, the Rl side chain lies at the heart of a crystal 
contact and is involved in numerous van-der-Waals 
(vdW) interactions and a hydrogen bond between the 
nitroxide oxygen and an ordered water molecule. This 
water molecule is again interacting with A119 of a 
neighboring azurin molecule (Figure 2B, D). For this Rl 
conformation (from here on referred to as Rl-II/IV) 
PISA calculates a solvent accessible surface area of 
260 A 2 , 171 A 2 of which are buried. The higher number 
of stabilizing contacts in Rl-II/IV is also reflected in a 
slightly higher quality of the electron density for this 
conformation (Figure 2 A, B). Interestingly, a glycerol 
molecule from the cryo-protectant interacts with the main 
chain atoms of Rl-II/IV (Additional file 1: Figure SI). In 
our structure, this molecule does not interact with 
the Rl side chain atoms. However, high concentrations of 



Mutual influence of Rl side chain and protein 
environment 

Figure 3A shows polar plots of the dihedral angles 
(xl- X^) measured from the two Rl side chains in 
Rl-I/III conformation. The two side chains were inde- 
pendently refined but have almost identical dihedral 
angles. The angles fit nicely to the angular distribution of 
a rotamer library derived from MD simulations of the free 
label [18]. In monomers II and IV, the crystal-packing 
environment precludes the Rl-I/III conformation due to 
steric constraints (Figure 4). The label adjusts to this new 
environment and therefore adopts the Rl-II/IV conform- 
ation. Again, the side chains were independently refined 
but show almost identical dihedral angles (Figure 3B). As 
described above, the Rl-II/IV conformation is stabilised 
by multiple vdW interactions and a bridging water mol- 
ecule (Figure 2D). Interestingly, the nearby K128 side 
chain rotates from the position found in monomers I and 
III by 180° to accommodate the Rl-II/IV conformation 
(Figure 4). Figure 3B shows polar plots of the dihedral 
angles measured in the Rl-II/IV conformation. Compared 
to Rl-I/III, the X2 and \ 4 angles change significantly. At 
the same time, Xs undergoes a smaller transition. Com- 
parison with the rotamer library reveals that %2> X4 an d to 
a lesser extend Xs of Rl-II/IV adjust to values that are less 
favorable for the free label (Figure 3B). 

Discussion 

In silico spin labelling programs are important tools to 
translate EPR-derived distances into structural informa- 
tion, despite rather unsatisfying deviations between ex- 
perimental and predicted distances [20,22,23]. To improve 
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Figure 2 Two different conformations of the R1 side chain. A) The R1 -l/lll conformation. The structure of monomer I is shown as green sticks; 
the neighbouring azurin monomer is coloured blue. The black mesh is the refined 2mFo-DFc electron density contoured at 1 o. The purple mesh 
is the positive difference electron density (mFo-DFc, contoured at 3.0 o) that was observed before the spin label was added to the structural 
model during refinement. B) The R1-II/IV conformation. The figure is analog to panel A, but in this case the structure of monomer II is shown in 
green. C) Interaction topology diagram (Ligplot+, [29]) of the R1 side chain in monomer I (R1 -l/lll conformation) and its environment. Covalent 
bonds are indicated by solid lines, polar interactions by dashed lines and vdW interactions by red arcs. D) Same as panel C but for the R1 side 
chain in monomer II (R1-II/IV conformation). 



the algorithms, it is important to understand the reasons 
for these deviations, and to identify the most promising 
points for improvements. In published benchmark studies, 
the performance of individual programs is usually com- 
pared based on a comparison of predicted and experimen- 
tal distances. But, the primary results of in silico spin 
labelling programs are predicted ensembles of spin labels. 
The two well-defined conformations of MTSSL in our 
crystal structure give rise to a number of inter-spin dis- 
tances in the crystal packing and we thought this to be an 
interesting way to investigate how well predicted ensem- 
bles correlate with predicted distances. It should be noted 
that spin label conformations found in a crystal structure 
can be biased by interactions with the crystalline protein 
environment. However, in contrast to crystals of small 
molecules, protein crystals are interspersed by large solv- 
ent channels and usually consist of around 50% solvent. 
Also, before X-ray data are collected at low temperature, 
the crystals are typically cryo-protected by soaking them 
in e.g. 35% glycerol (even higher concentrations are used 



for cryo-protection in PELDOR samples). This prevents 
the formation of crystalline ice in the solvent channels, 
which would otherwise destroy the crystal or severly de- 
grade its diffraction quality [30]. Thus, a spin label, which 
points into a solvent channel of a protein crystal is 
surrounded by the protein lattice, some ordered solvent 
molecules interacting with the spin-labelled protein and, a 
glassy, frozen solution of the solvent. Nevertheless, be- 
cause of the flash-cooling process and interactions with 
the protein lattice, the spin label conformations observed 
in a crystal might be different from those that are found at 
room temperature and/or in liquid solution. 

To address the issue of interactions between label and 
protein lattice in our analysis, the complete subset of the 
crystal packing shown in Figure 5A was selected as input 
model for two programs: MMM (rotamer approach) and 
mtsslWizard (accessible volume approach). This selection 
contains the spin labelled monomers I -IV (Figure 5 A) 
and, to complete the protein environment of each spin 
label, a symmetry equivalent of monomer IV (IV) was 
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Figure 3 Dihedral angles of the observed R1 side chains in azurin T21R1. A polar plot is shown for each x-angle. The angles measured in 
the R1 -l/lll conformation are plotted in A (coloured lines). The angles measured in the R1-II/IV conformation are plotted in B (coloured lines). The 
calculated angular distribution from the MTSSL rotamer library was digitized from [18] and is shown in each diagram (grey lines). 



additionally included (solvent molecules such as water 
and glycerol (see above) were excluded, because neither of 
the programs uses them). In principle, this selection can 
be seen as an experimental structure of one larger, conti- 
nous protein chain (termed Azurin B M> as in benchmark) 




R1-II/IV 



Figure 4 The protein environment forces the R1 side chain into 
the R1-II/IV conformation. A crystal contact of monomer II is 
shown in cyan. Side chains are represented as sticks. The alternative 
R1 -l/lll conformation of the R1 side-chain and the corresponding 
conformation of K128 are superposed (magenta). 



with multiple attached labels in either the Rl-I/III- or Rl- 
Il/IV-conformation. Importantly, the complete protein en- 
vironment of each label is present in this structure. For 
the purpose of the analysis below and for the reasons dis- 
cussed above, we assume that in a flash frozen PELDOR 
sample of Azurin BM , the spin label will adopt very similar 
conformations to those that we found in the flash cooled 
azurin T21R1 crystals (Figure 2). 

The table in Figure 5B compares the inter-label 
distances taken from this experimental structure to the 
equivalent distances obtained from mtsslWizard and 
MMM. Details of the MMM analysis (partition function, 
number of rotamers) are shown in Additional file 2: 
Figure S2. Further, a comparison of the predicted mean 
spin label positions is shown in Additional file 3: Figure 
S3. As found in larger distance-based benchmarks 
[20,22,23], both programs predict some of the distances 
quite accurately (e.g. II<->III), whereas large deviations 
are found for other distances (e.g. III<-^IV'). This indi- 
cates that errors stemming from the generation of the 
spin label ensembles are sometimes compensated for by 
the relative geometric arrangement of a pair of ensem- 
bles. A similar observation has also been made in the 
crystal structure of the Spal5 chaperone [31]. Figure 5C 
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1 


II 


III 


IV 


1 




77.9 
73.0 
77.0 


56.0 
57.9 
60.8 


61.2 
67.0 
65.3 


II 






33.0 
34.9 
33.2 


25.6 
21.9 
21.3 


III 






. 


7.7 
16.2 
12.4 


IV 











R-II/IV 
monomer II 



R-l/lll 
monomer I 





Figure 5 Experimental and predicted distances between different azurin monomers in the azurin T21R1 crystal structure. A) Selection 
from the crystal packing that was used as input for the spin labelling programs. Red spheres indicate the label positions. The grey monomer was 
included to complete the environment of the spin label on monomer I. B) Distances between the monomers in panel A). Distances from the 
X-ray structure are pink, distances predicted by MMM are green and distances predicted by mtsslWizard are blue. For MMM the distances were 
taken from the .html output file generated by the program. C) close up of the ll-lll distance for mtsslWizard (top) and MMM (bottom). The 
experimental distance was taken from the X-ray structure and is shown in pink. The distance vector predicted by mtsslWizard is shown in blue 
and connects the geometric averages of the two ensembles predicted by the program (blue sticks). The green vector is the distance vector 
between the occupancy weighted geometric average of the MMM predicted ensembles (green sticks). For MMM, the occupancy of the individual 
rotamers is represented by the thickness of the sticks. Thus, only rotamers with a predicted occupancy larger than zero are shown. 



illustrates this for the II-III distance: Whereas the abso- 
lute values of the distance vectors are very similar 
between X-ray structure, MMM and mtsslWizard, their 
directions differ considerably. This pair of labels is ana- 
lysed in more detail below. 

The Rl-I/III label in Figure 5C is located at a relatively 
"open" labelling site on monomer III (Figure 5A). In the 
crystal structure, the spin label at this site adopts a 
conformation that fits rather well to calculated rotamer 
libraries (Figure 3A). Figure 6 A (right) shows a close-up 
of this site with the MTSSL ensemble predicted by 
MMM. Indeed, it contains some rotamers that are simi- 
lar to the conformation found in the X-ray structure. 
But, even when only the rotamers from this ensemble 
are considered, the nitroxide group can still sample a 



large volume (Figure 6A, right). Thus, especially at open 
sites, it is crucial to also correctly predict the occupancy 
of the individual rotamers. In Figure 6A occupancies 
predicted by MMM are illustrated by the radius of the 
stick models [18]. In this case, the program outputs very 
similar occupancies for all rotamers. In contrast, the X-ray 
structure exhibits only one defined conformation and no 
indications for further conformations were found in the 
electron density (Figure 2A). Assuming, that due to its 
interactions with the environment (Figure 2A), this single 
conformation of the spin label would also be present in a 
PELDOR sample of Azurin BM , the failure to correctly pre- 
dict the occupancy of the rotamers effectively neutralizes 
the advantage of the rotamer approach. This becomes 
clear by comparing the MMM ensemble with an ensemble 
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x-ray 

mtsslWizard 



Figure 6 The azurin T21R1 X-ray structure (MTSSL: pink sticks) in comparison with MTSSL rotamers predicted by MMM (A, green) and 
mtsslWizard (B, blue). The predicted ensembles are shown as sticks and for MMM (A) the predicted occupancy is represented by the thickness 
of the sticks. Thus, only rotamers with a predicted occupancy larger than zero are shown. 



predicted by mtsslWizard at the same site (Figure 6B, 
right). The overall volumes of both predicted ensembles 
are similar. Thereby, the uncertainty that will be intro- 
duced into predicted distances from or to this site will be 
similar between rotamer and accessible volume approach. 

Clearly, the influence of the protein environment will 
be more pronounced at tight labelling sites, such as the 
one on monomer II (Figure 5C, left). As a result, devia- 
tions from the rotamer libraries become more likely, and 
we observe this in the crystal structure (Figure 3B). Con- 
sequentially, the Rl-II/IV conformation cannot be cor- 
rectly predicted when the rotamer library shown in 
Figure 3 is used (Figure 6 A, left). It should be noted that 
this limitation of the rotamer approach at tight sites was 
pointed out by the authors of the MMM software [32]. 
In such cases, the rotamer approach effectively boils 
down to an accessible volume approach and would again 
deliver results that are very similar to the latter if the 
occupancies for all rotamers were set to the same value 
(Figure 6, left). 

In the examples above, the spin labelled X-ray struc- 
ture was used as a basis for in silico spin labelling. 
Usually, the problem is even more difficult, since it is 
unknown, how the protein will structurally react to the 



addition of the spin label. The K128 side chain in our 
structure is an example for a structural response of the 
protein (Figure 4). In silico spin labelling programs try to 
deal with this problem by allowing a certain number of 
clashes between protein and label. However, by doing 
this, the ensemble of created rotamers will simply grow 
in size whereas in reality the ensemble might merely 
change its shape, not necessarily its size. In the end this 
again leads to an increased uncertainty of the prediction. 

Conclusions 

Our observations vividly illustrate why in our test case 
(and possibly also in general), the accessible volume 
approach and the more sophisticated rotamer approach 
often deliver very similar results. In essence, employing 
rotamer libraries will only increase the accuracy, when 
not only the rotamers but also their occupancy can be 
correctly predicted. The occupancy prediction is cru- 
cially dependent on the interaction of the label with its 
environment. The current software programs use only 
relatively simple (but fast) descriptions of the protein 
environment, whereas solvent molecules such as the 
bridging water shown in Figure 2 are completely ignored. 
Recently, the formation of hydrogen bonds has been 
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found to be a very common type of interaction for nitrox- 
ide spin labels on proteins [33]. Also ignored are cryo 
protectants such as ethylene glycole or glycerol which are 
used at high percentages in PELDOR samples. These can 
bind close to the labelling site and thereby potentially 
influence the label dynamics (see above and Additional 
file 1: Figure SI). For labelling sites close to lipid bilayers 
or detergent micelles it is also important to consider label- 
lipid or label-detergent interactions since these will likely 
have profound effects on the conformation of the label. 
Thus, to increase the prediction accuracy, more sophisti- 
cated algorithms, which accurately account for label- 
environment interactions will have to be employed. It 
has been tried multiple times to use MD simulations 
for this purpose, but so far, the increased effort does 
not seem to pay off in terms of better prediction ac- 
curacy [34]. Promising ways to alleviate the described 
difficulties on the experimental side would be the use 
of spin labels with shorter or conformationally restrained 
linkers, such as the RX side chain [35] or the recently pub- 
lished VI side chain [36]. 

Methods 

Cloning, protein expression, purification and spin 
labelling 

The gene for azurin (azu, PA4922) was PCR amplified 
from genomic Pseudomonas aeruginosa DNA using the 
PCR primers 5 -TTATAACCATGGCCGAGTGCTCGG 
TGG-3' and 5 -CACCCTGACCCTGAAGTGAGAGCT 
CTTATAA-3'. The resulting PCR product did not con- 
tain the coding region for the N-terminal signal peptide 
of azurin (residues - 20 - 0), so that the target protein 
(residues 1 - 128 of azurin) could be expressed intracel- 
lularly in E. coli. The PCR product was then cloned into 
the vector pEHISGSTTEV (Huanting Liu, Biomedical 
Sciences Research Center, University of St Andrews, 
UK) via restriction enzymes Ncol and SacI, resulting in 
an expression construct with a TEV cleavable N-terminal 
His6-GST (glutathione S -transferase) tag. The T21C 
mutant was introduced into this construct using PCR. 
The resulting construct was transformed into E. coli 
Rosetta cells. A single colony was picked and grown over 
night in 50 ml of 2xYeast-Trypton (2YT) media supple- 
mented with 100 (ig/ml ampicillin and 17 (ig/ml chlor- 
amphenicol with shaking at 37°C. On the following 
day, 1 1 of 2YT medium with 50 (ig/ml kanamycin and 
17 (ig/ml chloramphenicol were inoculated with 20 ml of 
the overnight culture and grown to an OD600 of 1.0. 
Protein expression was then induced by addition of 
0.3 mM iso-propoyl-beta-thiogalactoside (IPTG). The 
protein expression was allowed to proceed for 3 h at 
37°C with shaking at 200 rpm. The cells were then 
harvested by centrifugation at 2800 g, resuspended in 
100 ml of lysis buffer (20 mM Tris-HCl pH 7.5, 500 mM 



NaCl, 30 mM Imidazol) and lysed with a cell disrupter at 
30 kPsi (Constant Systems). Cell debris and insoluble pro- 
teins were spun down at 32.000 g for 15 min at 4°C. The 
soluble fraction was mixed with 1.5 ml Ni-NTA resin 
(Quiagen, pre-equilibrated in lysis buffer) and incubated 
for 1 h at 4°C with shaking. The resin was washed with 
100 ml of lysis buffer, followed by 50 ml of lysis buffer 
supplemented with 1 mM tris(2-carboxyethyl)phosphine 
(TCEP) to reduce the introduced cysteine residue. The re- 
ducing agent was then quickly removed by washing the 
column with 50 ml of lysis buffer, directly followed by 
addition of 15 ml elution buffer (20 mM Tris-HCl pH 7.5, 
500 mM NaCl, 1 M Imidazol) containing 0.7 mM of 
MTSSL. A large excess (~20x) of MTSSL was used, since 
the GST-tag of the expression construct also contained 
four cysteine residues. The labeling reaction was trans- 
ferred to dialysis tubing and dialyzed over night against 5 1 
of dialysis buffer (20 mM Tris-HCl pH 7.5, 500 mM 
NaCl). On the next day, 4 mg TEV protease were added 
to the sample to cleave the GST- tag. The cleavage reaction 
was incubated for 3 h at room temperature. The sample 
was then concentrated to a volume of 2 ml, supplemented 
with 1 mM CuCl2 and loaded onto a Superdex200 
16/60 column (GE) equilibrated with gel filtration buffer 
(10 mM Tris-HCl pH 8.0, 150 mM NaCl). Labelled mono- 
meric azurin eluted at a volume of -100 ml and had an 
intense blue colour. 



Crystallisation, data collection and refinement 

Purified and MTSSL-labelled azurin T21R1 was concen- 
trated to 10 mg/ml for crystallization and sitting drop 
crystallization setups were prepared with the commercial 
JCSG + screen (Molecular Dimensions) in MRC plates (Mo- 
lecular Dimensions). Blue, plate shaped azurin crystals grew 
within 2-3 days at room temperature in conditions A2, B7, 
C4 and Cll. The components of these conditions where 
then used for a stochastic optimization and led to the final 
condition: 2.14 M ammonium sulfate, 0.28 M ammonium 
nitrate, 0.1 M sodium cacodylate pH 6.0. For data collection, 
the crystals were harvested and cryo-protected with 35% 
glycerol prior to flash cooling in liquid nitrogen. A diffrac- 
tion data set was collected at BESSYII, BL14.1 (Berlin) using 
a PILATUS 6M detector. The data were indexed 
(iMOSFLM [37]) in space group PI and processed with 
iMOSFLM, POINTLESS and SCALA [38]. The structure of 
Zn 2+ -bound azurin (PDB-ID: 1E67) was used as model for 
molecular replacement with PHASER [39]. The program lo- 
cated all four azurin molecules in the asymmetric unit. The 
monomers (I-IV) are related by a two-fold non-crystallo- 
graphic symmetry between I, II and III, IV. The model 
was refined automatically using PHENIX.REFINE [40] 
and by hand using COOT [41]. Data collection and refine- 
ment statistics are listed in Table 1. 
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Additional files 



Additional file 1: Figure SI. A glycerol molecule (lilac) from the cryo 
protectant interacts (dotted lines) with the main chain atoms of the R1 
label. The protein structure is shown as green sticks, the neighbouring 
protein in the crystal is shown as blue sticks. Water molecules are shown 
as red crosses. 

Additional file 2: Figure S2. Cartoon image of the structure that was 
used as input for the in silico spin labelling programs. The colouring is 
identical to Figure 5 but the partition functions (p.f.) and number of 
rotamers (#) generated by MMM are given. The label on monomer IV was 
ignored because its complete protein environment is not part of the 
ensemble. 

Additional file 3: Figure S3. Cartoon image of the structure that was 
used as input for the in silico spin labelling programs. The colouring is 
identical to Figure 5. The average coordinates from MMM (green spheres), 
mtsslWizard (blue spheres) and the N1 coordinates from the X-ray structure 
are shown. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

NF & GH performed research. GH and OS designed research. GH and OS 
wrote the paper. All authors read and approved the final manuscript. 

Acknowledgements 

We thank Drs Uwe Muller and Manfred Weiss for granting us beam time 
on BL14.1, BESSY II, Berlin. We would like to thank Dr Hartmut Niemann for 
access to his X-ray source. Funding from DFG (SFB 813, project A6) is 
gratefully acknowledged. The pEHisGS^EV vector was a gift from Huanting 
Liu from the University of St Andrews, UK. 

Received: 7 March 2014 Accepted: 8 May 2014 
Published: 29 May 2014 

References 

1. Reginsson GW, Schiemann 0: Studying biomolecular complexes with 
pulsed electron-electron double resonance spectroscopy. Biochem Soc 
Trans 2011,39:128-139. 

2. Jeschke G: DEER distance measurements on proteins. Annu Rev Phys Chem 
2012, 63:419-446. 

3. Schiemann 0, Prisner TF: Long-range distance determinations in 
biomacromolecules by EPR spectroscopy. Q Rev Biophys 2007, 40:1-53. 

4. Berliner LJ, Grunwald J, Hankovszky HO, Hideg K: A novel reversible 
thiol-specific spin label: papain active site labeling and inhibition. 

Anal Biochem 1982, 119:450-455. 

5. Klare JP, Steinhoff HJ: Spin labeling EPR. Photosynth Res 2009, 102:377-390. 

6. Altenbach C, Marti T, Khorana HG, Hubbell WL: Transmembrane protein 
structure: spin labeling of bacteriorhodopsin mutants. Science 1990, 
248:1088-1092. 

7. Hubbell WL, Lopez CJ, Altenbach C, Yang Z: Technological advances in 
site-directed spin labeling of proteins. Curr Opin Struct Biol 2013, 
23(5)725-733. doi:1 0.1 01 6/j.sbi.201 3.06.008. 

8. Smirnova I, Kasho V, Choe J-Y, Altenbach C, Hubbell WL, Kaback HR: Sugar 
binding induces an outward facing conformation of LacY. Proc Natl Acad 
SciUSA 2007, 104:16504-16509. 

9. Altenbach C, Kusnetzow AK, Ernst OP, Hofmann KP, Hubbell WL: 
High-resolution distance mapping in rhodopsin reveals the pattern of 
helix movement due to activation. Proc Natl Acad Sci USA 2008, 
105:7439-7444. 

10. Ward R, Pliotas C, Branigan E, Hacker C, Rasmussen A, Hagelueken G, 
Booth IR, Miller S, Lucocq J, Naismith JH, Schiemann O: Probing the 
structure of the mechanosensitive channel of small conductance in lipid 
bilayers with pulsed electron-electron double resonance. Biophys J 2014, 
106:834-842. 

11. Hilger D, Polyhach Y, Padan E, Jung H, Jeschke G: High-resolution structure 
of a Na+/H + antiporter dimer obtained by pulsed electron 



paramagnetic resonance distance measurements. Biophys J 2007, 
93:3675-3683. 

12. Park S-Y, Borbat PP, Gonzalez-Bonet G, Bhatnagar J, Pollard AM, Freed JH, 
Bilwes AM, Crane BR: Reconstruction of the chemotaxis receptor-kinase 
assembly. Nat Struct Mol Biol 2006, 13:400-407. 

13. Langen R, Oh KJ, Cascio D, Hubbell WL: Crystal structures of spin labeled 
T4 lysozyme mutants: implications for the interpretation of EPR spectra 
in terms of structure. Biochemistry 2000, 39:8396-8405. 

14. Guo Z, Cascio D, Hideg K, Hubbell W: Structural determinants of nitroxide 
motion in spin-labeled proteins: solvent-exposed sites in helix B of T4 
lysozyme. Protein Sci 2008, 1 7:228-239. 

15. Guo Z, Cascio D, Hideg K, Kalai T, Hubbell WL: Structural determinants 
of nitroxide motion in spin-labeled proteins: tertiary contact and 
solvent-inaccessible sites in helix G of T4 lysozyme. Protein Sci 2007, 
16:1069-1086. 

16. Fleissner MR, Cascio D, Hubbell WL: Structural origin of weakly ordered 
nitroxide motion in spin-labeled proteins. Protein Sci 2009, 18:893-908. 

17. Sale K, Song L, Liu Y-S, Perozo E, Fajer PG: Explicit treatment of spin labels 
in modeling of distance constraints from dipolar EPR and DEER. 

J Am Chem Soc 2005, 127:9334-9335. 

18. Polyhach Y, Bordignon E, Jeschke G: Rotamer libraries of spin labelled 
cysteines for protein studies. Phys Chem Chem Phys 201 1, 13:2356-2366. 

19. Hatmal MM, Li Y, Hegde BG, Hegde PB, Jao CC, Langen R, Haworth IS: 
Computer modeling of nitroxide spin labels on proteins. Biopolymers 
2011,97:35-44. 

20. Hagelueken G, Ward R, Naismith JH, Schiemann O: MtsslWizard: in silico 
spin-labeling and generation of distance distributions in PyMOL. 
Appl Magn Reson 2012, 42:377-391. 

21. Hagelueken G, Abdullin D, Ward R: mtsslSuite: In silico spin labelling, 
trilateration and distance-constrained rigid body docking in PyMOL. 
Mol Phys 2013, 111:2757-2766. doi:1 0.1 080/00268976.201 3.809804. 

22. Alexander NS, Stein RA, Koteiche HA, Kaufmann KW, Mchaourab HS, Meiler 
J: RosettaEPR: rotamer library for spin label structure and dynamics. 
PLoS One 2013, 8:e72851. 

23. Jeschke G: Conformational dynamics and distribution of nitroxide spin 
labels. Prog Nucl Magn Reson Spectrosc 2013, 72:42-60. 

24. Karlsson BG, PascherT, Nordling M, Arvidsson RH, Lundberg LG: Expression 
of the blue copper protein azurin from Pseudomonas aeruginosa in 
Escherichia coli. FEBS Lett 1989, 246:21 1-217. 

25. Liu H, Naismith JH: An efficient one-step site-directed deletion, insertion, 
single and multiple-site plasmid mutagenesis protocol. BMC Biotechnol 
2008, 8:91. 

26. Nar H, Huber R, Messerschmidt A, Filippou AC, Barth M, Jaquinod M, van de 
Kamp M, Canters GW: Characterization and crystal structure of zinc azurin, a 
by-product of heterologous expression in Escherichia coli of Pseudomonas 
aeruginosa copper azurin. Eur J Biochem 1992, 205:1 123-1 129. 

27. Sauter NK, Grosse-Kunstleve RW, Adams PD: Improved statistics for 
determining the Patterson symmetry from unmerged diffraction 
intensities. J Appl Crystallogr 2006, 39:158-168. 

28. Krissinel E, Henrick K: Inference of macromolecular assemblies from 
crystalline state. J Mol Biol 2007, 372:774-797. 

29. Laskowski RA, Swindells MB: LigPlot+: multiple ligand-protein 
interaction diagrams for drug discovery. J Chem Inf Model 2011, 
51:2778-2786. 

30. Garman E: Cool data: quantity AND quality. Acta Crystallogr D Biol 
Crystallogr 1999, 55:1641-1653. 

31. Lillington JED, Lovett JE, Johnson S, Roversi P, Timmel CR, Lea SM: Shigella 
flexneri Spa15 crystal structure verified in solution by double electron 
electron resonance. J Mol Biol 201 1, 405:427-435. 

32. Polyhach Y, Jeschke G: Prediction of favourable sites for spin labelling of 
proteins. J Spectroscopy 2010, 24:651-659. 

33. Urban L, Steinhoff HJ: Hydrogen bonding to the nitroxide of protein 
bound spin labels. Mol Phys 2013, 1 1 1:2873-2881. doi:1 0.1 080/ 
00268976.2013.804217. 

34. Klose D, Klare JP, Grohmann D, Kay CWM, Werner F, Steinhoff HJ: 
Simulation vs. reality: a comparison of in silico distance predictions with 
DEER and FRET measurements. PLoS One 2012, 7:e39492. 

35. Fleissner MR, Bridges MD, Brooks EK, Cascio D, Kalai T, Hideg K, Hubbell WL: 
Structure and dynamics of a conformationally constrained nitroxide side 
chain and applications in EPR spectroscopy. Proc Natl Acad Sci USA 201 1, 
108:16241-16246. 



Florin et al. BMC Structural Biology 2014, 14:16 
http://www.biomedcentral.eom/1 472-6807/1 4/1 6 



Page 10 of 10 



36. Toledo Warshaviak D, Khramtsov W, Cascio D, Altenbach C, Hubbell WL: 
Structure and dynamics of an imidazoline nitroxide side chain with 
strongly hindered internal motion in proteins. J Magn Reson 2013, 

232:53-61. 

37. Battye TGG, Kontogiannis L, Johnson 0, Powell HR, Leslie AGW: iMOSFLM: 

a new graphical interface for diffraction-image processing with MOSFLM. 

Acta Crystallogr D Biol Crystallogr 201 1 , 67:271 -281 . 

38. Evans P: Scaling and assessment of data quality. Acta Crystallogr D Biol 
Crystallogr 2005, 62:72-82. 

39. McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, Read RJ: 
Phasercrystallographic software. J Appl Crystallogr 2007, 40:658-674. 

40. Adams P, Grosse-Kunstleve R, Hung L, loerger T, McCoy A, Moriarty N, Read R, 
Sacchettini J, Sauter N, Terwilliger T: PHENIX: building new software for 
automated crystallographic structure determination. Acta Crystallogr D Biol 
Crystallogr 2002, 58:1 948-1 954. 

41. Emsley P, Cowtan K: Coot: model-building tools for molecular graphics. 
Acta Crystallogr D Biol Crystallogr 2004, 60:21 26-21 32. 



doi:1 0.1 1 86/1 472-6807-1 4-1 6 

Cite this article as: Florin et al.: High-resolution crystal structure of spin 
labelled (T21R1) azurin from Pseudomonas aeruginosa: a challenging 
structural benchmark for in silico spin labelling algorithms. BMC 

Structural Biology 2014 14:16. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f~\ RiftMM i rpntral 

www.biomedcentral.com/submit momea central 



